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ABSTRACT 

Plutinos are Kuiper-belt objects that share the 3:2 Neptune resonance with 
Pluto. The long-term stability of Plutino orbits depends on their eccentric- 
ity. Plutinos with eccentricities close to Pluto (fractional eccentricity difference 
Ae/e p = \e — e p \/e p & 0.1) can be stable because the longitude difference librates, 
in a manner similar to the tadpole and horseshoe libration in coorbital satellites. 
Plutinos with Ae/e p < 0.3 can also be stable; the longitude difference circulates 
and close encounters are possible, but the effects of Pluto are weak because the 
encounter velocity is high. Orbits with intermediate eccentricity differences are 
likely to be unstable over the age of the solar system, in the sense that encoun- 
ters with Pluto drive them out of the 3:2 Neptune resonance and thus into close 
encounters with Neptune. This mechanism may be a source of Jupiter-family 
comets. 



Subject headings: planets and satellites: Pluto — Kuiper Belt, Oort cloud 
celestial mechanics, stellar dynamics 



Introduction 



The orbit of Pluto has a number of unusual features. It has the highest eccentricity 
(e p = 0.253) and inclination (i p = 17.1°) of any planet in the solar system. It crosses 
Neptune's orbit and hence is susceptible to strong perturbations during close encounters 
with that planet. However, close encounters do not occur because Pluto is locked into a 
3:2 orbital resonance with Neptune, which ensures that conjunctions occur near Pluto's 
aphelion ( pohen fc Hubbard 1965| ). More precisely, the critical argument 3A P — 2A n — w v 
librates around 180° with a period of 1.99 x 10 4 yr and an amplitude of 82°; here \ p and A n are 
the mean longitudes of Pluto and Neptune and w p is Pluto's longitude of perihelion. Other 
resonances are present with longer periods: for example, Pluto's argument of perihelion 
librates with an amplitude of 23° and a period of 3.8 Myr. In part because of its rich set 
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of resonances, Pluto's orbit is chaotic, although it exhibits no large-scale irregular behavior 
over Gyr timescales (see |Malhotra fc Williams 1997| for a comprehensive review of Pluto's 



orbit). For reference, Pluto's semimajor axis and orbital period are a p = 39.774 AU and 
P p = 250.85 yr. 

The most compelling explanation for Pluto's remarkable orbit was given by Malhotra 
(1993, 1995). Malhotra argues that Pluto formed in a low-eccentricity, low-inclination orbit in 
the protoplanetary disk beyond Neptune. Subsequent gravitational scattering and ejection of 
planetesimals in the disk by all four giant planets caused Neptune's orbit to migrate outwards 
( |Fcrnandez fc Ip 1984| ). As its orbit expands, Neptune's orbital resonances sweep through 
the disk, first capturing Pluto into the 3:2 resonance and then pumping up its eccentricity. If 
Pluto's orbit was circular before capture, its present eccentricity implies that it was captured 
when Neptune's semimajor axis was 0.814 times its current value or 24.6 AU (eq. |5|). This 
process may also excite Pluto's inclination although the details are less certain ( Malhotra 
1998ft . 



Malhotra's argument predicts that most Kuiper belt objects with 30 AU ^ a ^ 50 AU 
should also be captured — and presently located — in Neptune resonances ( Malhotra 1995| ). 
This prediction has proved to be correct: of the ~ 90 Kuiper belt objects with reliable 
orbits as of 1999 January 1, over 30% have semimajor axes within 1% of the 3:2 resonance 
(although this number is exaggerated by observational selection effects). These objects have 
come to be called Plutinos, since they share the 3:2 resonance with Pluto (see Malhotra et 



al. 1999| for a recent review of the Kuiper belt). 



Almost all studies of the dynamics of the Kuiper belt so far have neglected the gravi- 
tational influence of Pluto, because of its small mass (M p /M Q = 7.40 x 10~ 9 for the Pluto- 



Charon system, [Stern 19TJ2], |lholen fc Buie 1997[). However, like the Trojan asteroids and 



Jupiter, or the Saturn coorbital satellites Janus and Epimetheus, Pluto and the Plutinos 
share a common semimajor axis and hence even the weak gravitational force from Pluto 
can have a substantial influence on the longitude of a Plutino relative to Pluto. A crude 
illustration of the importance of Pluto's gravity is to note that the half-width of the 3:2 
resonance, (Aa/a) res ~ 0.01 for 0.2 ^ e ^ 0.3 (the maximum fractional amplitude of stable 
libration in semi- major axis, Malhotra 1996| ), is only a few times larger than the Hill radius 



of Pluto, (Aa/o) H = (M p /3M G y/ 3 = 0.0014. 

The goal of this paper is to explore the dynamical interactions between Pluto and 
Plutinos and their consequences for the present structure of the Kuiper belt. Section 2 
provides an approximate analytical description of the interactions, §3 describes the results 
of numerical orbit integrations, and §4 contains a discussion. 
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2. Analysis 

We examine a simplified model solar system containing only the Sun and Neptune, 
with masses M and M n ; we assume that Neptune's orbit is circular and neglect all orbital 
inclinations. We describe the motion of the Plutino using the canonical variables 

Xl = {GM Q a) l '\ Vl = A, 

x 2 = (GM Q a) 1/2 [l - (1 - e 2 ) 1/2 ], 2/2 = -w, (1) 

where a and e are the semi-major axis and eccentricity, A and w are the mean longitude and 
longitude of perihelion. The same variables for Neptune or Pluto are denoted by adding a 
subscript "n" or "p". 

We consider only the resonant perturbations exerted by Neptune, which can only depend 
on angles as an even function of the combination 3|/i — 1y n \ + U2- Thus the Hamiltonian of 
a Plutino may be written 

if (x, y, t) = H K { Xl ) + A(x, 3j/i - 2y nl + y 2 ) + £(x, y, t), (2) 

where Hk(xi) = — \{GM & ) 2 /x\ is the Kepler Hamiltonian, A is the resonant potential from 
Neptune, and _£>(x, y,t) is the potential from Pluto. The same Hamiltonian describes the 
motion of Pluto if we set B = 0. 

Now impose a canonical transformation to new variables (J, w) defined by the generating 
function 

^(Jjy,*) = Ji(3j/i - 2y nl + y 2 ) + \J2hj2 -Vja)- (3) 

Thus wi = dS/dJi = 3y 1 - 2y nl +y 2 , w 2 = dS/dJ 2 = \{y 2 -y p2 ), while x x = dS/dy 1 = 3Ji, 
x 2 = dS/dy 2 — J± + I J 2 , and the new Hamiltonian is 

dS 

H(J, w ,t) = H + — = ^(3Ji) - 2y nl Jx - iy p2 J 2 + A(3Ji, Ji + \ J 2 , ^i) + 5, (4) 

where y„i is the mean motion of Neptune and — y p2 is the apsidal precession rate of Pluto. 

To proceed further we make use of the fact that B/A = 0(M p /M n ) <^ 1. Thus we can 
divide the motion of a Plutino into "fast" and "slow" parts. The fast motion is determined 
by the Kepler Hamiltonian Hk and the resonant potential A (this is the opposite of normal 
usage, where the resonant perturbations from Neptune are regarded as "slow" compared to 
non-resonant perturbations). The slow variations are caused by the Pluto potential B. 

First we examine the fast motion. We drop the potential B, so that w 2 is ignorable and 
J 2 = 3x 2 - X! = (GM Q a) 1/2 [2 - 3(1 - e 2 ) 1 / 2 ] is a constant of the motion. Thus if Pluto was 
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captured into resonance from a circular orbit with semimajor axis a^, its present semimajor 
axis and eccentricity are related by 

e 2 = I - Ka./a) 1 / 2 - i(a,/a). (5) 

We write J\ = J\ r + A J\ where J\ r is chosen to satisfy the resonance condition for the 
Kepler Hamiltonian, 

2y nl = 3 yi = 3— or J lr = ^ = p-^. (6) 

Since A/H K = O(M n /M ) < 1 we expect |AJi| <C Ji r , so we can expand H K to second 
order in AJi; dropping unimportant constant terms the fast motion is determined by the 
Hamiltonian 

HfiAJ^wx) = \{^§jr ) ) (AJ 1 ) 2 + A(3J lr + 3AJ 1 ,J lr + AJ l + iJ 2 ,w 1 ) 
27 

(AJ X ) 2 + A(3J lr + 3AJ 1 ,J lr + AJ 1 + U 2 ,w 1 ). (7) 



2a 

The Hamiltonian is autonomous and hence has a conserved energy Ef = Hf and action 
/ = (27r) _1 § AJidwi. The motion is along the level surfaces of Hf in the (AJi,u>i) plane 
and typically consists of either libration {w\ oscillates between fixed limits) or circulation {w\ 
increases or decreases without reversing), just as in the case of the pendulum Hamiltonian. 
The stable equilibrium solutions (i.e. zero-amplitude libration) are given by 



a 2 (dA 1 dA\ dA d 2 A 

9 V dxi 3 dx 2 1 ' dwi dwj 



The slow motion is determined by averaging the Hamiltonian (^) over the fast motion: 

H s (J 2 ,w 2 ,t) = E f {J 2 ) - \y p2 J 2 + (B). (9) 

Here (■) indicates an average over one period of the fast motion. The fast energy Ef depends 
on J 2 through the constraint that the fast action / is adiabatically invariant. 



2.1. Solutions with zero-amplitude libration 



The solutions to the fast and slow equations of motion are particularly simple in the case 
where the fast libration amplitude is zero for both Pluto and the Plutino. This approximation 
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is not particularly realistic — the libration amplitude of Pluto is 82° — but illustrates the 
principal features of the Plutino motions. 

In this case the fast energy is 

E f = A(3J lr ,J lr + %J 2 ,w lr ) (10) 

where w Xr is the equilibrium angle given by equation (j8|) and we have dropped much smaller 
terms that are 0(A 2 ). For simplicity we shall assume that there is only one stable equilibrium 
point, that is, one solution to equations (§) for given J 2 . At the equilibrium point the fast 
action is / = 0, and the slow Hamiltonian (|9|) is 

H s { J 2 , w 2 , t) = A(3J lr , J lr + § J 2 , w lr ) - \y p2 J 2 + (B(3J ln J Xr + § J 2 , y x , y 2 , t)) + 0{A 2 ), (11) 

where y\ = \{wi r + 2y n i — y 2 ). Since Pluto also is assumed to have zero libration amplitude, 
y P i = l(w lr + 2y nl - y p2 ). Thus 

3(y P i - yi) = y% - y P 2 = 3w 2 mod(2?r); (12) 

that is, the difference in longitude of perihelion between Pluto and the Plutino is three times 
the difference in mean longitude. The same result will hold true on average even if the 
libration amplitudes are non-zero. 

Now let us assume in addition that the eccentricities of both Pluto and the Plutino are 
small. Since their semi-major axes are the same, the gravitational potential from Pluto at 
the Plutino may be written 



B = GM P 



r r p 1 \ _ GM P 
Jr p | 3 |r-r p |y a 



cos(yi -y pX ) 



2 1 sin 5 (yi -y P i) 



(13) 



Using equation (]T2|) this simplifies to 

u ( 

2\ sin ^2 



B = 9^( mwt . _2_ V (14) 



Moreover 

% = |^(3J ln J lr + U p2 , w lr ) + 0(A 2 ) 



dt dx 2 ^" ±r, " XT ' 3 p2 

thus the slow Hamiltonian can be rewritten as 

dA 
'dx 2 



(15) 



H S (J 2 , W 2 ) = A(3J lr , J lr + | J 2 , Wlr) - |^2^ (3 Ji r , Jl r + | J p 2, W X 

ox 2 

GM r , ( 1 



f r^£ I COSW2 _ - ] + (A 2 ). (16) 
a \ 2 1 sin ±w 2 \ 
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The interesting behaviour occurs when the actions J 2 and J p2 are similar (i.e. the eccen- 
tricities of Pluto and the Plutino are similar), so we write J 2 = J p2 + AJ 2 and expand A to 
second order in AJ 2 . Dropping unimportant constants and terms of 0(A 2 ) we have 



H s (AJ 2 ,w 2 ) = ^AJ 2 v4 22 (3Ji r , J lr + \ J P 2,w lr ) + ^ 



COS W2 - . . \ r J . (17) 

2 1 sin~w 2 |/ 



where A 22 = <9 2 t4/<9x 2 . 



This Hamiltonian is strongly reminiscent of the Hamiltonian for a test particle coorbiting 
with a satellite, 

.wax 3 2 GM ( 1 \ 

H c (Ax,w) = --^Ax 2 + [cosw- . 1 ; 18 

2cr a y 2|sm|io|/ 

here Ax = x\ — x s i and w = \\ — X s i are conjugate variables, and the subscript s denotes 
orbital elements of the satellite. In this case the torques from the satellite lead to changes in 
semi-major axis; for a Plutino the semi-major axis is locked to Neptune's by the resonance, 
so torques from Pluto lead to changes in the eccentricity instead. 

When v4 22 < many of the features of orbits in the slow Hamiltonian (|17|) follow 
immediately from the analogy with the Hamiltonian fll8|), which has been studied by many 
authors (e.g. Yoder et al. 1983, Namouni et al. 1999). The trajectories are determined by 
the level surfaces of the Hamiltonian. The equilibrium solutions correspond to the triangular 
Lagrange points in the coorbital case: AJ 2 = 0, w 2 = ±60°, H s = — |GM p / 'a; the eccentricity 
of the Plutino equals the eccentricity of Pluto, the mean longitude leads or lags by 60°, and 
the perihelia are 180° apart. These solutions are maxima of the potential from Pluto. For 
smaller values of H s , the orbits librate around the triangular points ("tadpole orbits"). 
Small-amplitude tadpole librations have frequency u given by 

9 1 „ GM„ . . 

uj 2 = —A22 P ~. 19 

4 a 

The tadpole orbits merge at the separatrix orbit, H s = —^GM P / a; for this orbit the minimum 
separation is w 2 , m in = 23.91°. Even smaller values of H s yield "horseshoe" orbits, with 
turning points at ±w 2jm i n where H s = (GM p / a) (cos w 2 , m m — 1| sin |w 2jm i n | _1 ). For all tadpole 
and horseshoe orbits, the maximum and minimum values of AJ 2 occur at u> 2 = ±60°, and 
are given by 

18 / GM„^ 1/2 



AJ 2 = ± 



(20) 



.A 2 2 V 2a 

Eventually the theory breaks down, when u> 2jm i n is small enough that adiabatic invariance is 
no longer a valid approximation. 
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2.2. The resonant potential from Neptune 

For quantitative applications we must evaluate the resonant Neptune potential A(x, Wi). 
For small eccentricities, the potential can be derived analytically, 



A( \ GM « 
A(x,u;i) = 



\bf l2 {ct) + \e 2 (2aD + a 2 D 2 )bfJ 2 (a) + §e(5 + aD)bfj 2 {a) cosw x 



+ |e 2 (104 + 22aD + a 2 D 2 )bf/ 2 {a) cos 27^ + 0(e 3 



(21) 



where a = a n /a < 1, D = d/da, 

b U) (a) _ I f 2n #cosj0 

bl / 2[a) ~nJo (l-2ocos0 + o 2 ) 1 /2 ^ 

is a Laplace coefficient, Neptune is assumed to be on a circular orbit, and inclinations are 
neglected. 

Unfortunately, in the present case the high eccentricity of the Plutino orbits makes this 
expansion invalid. However, we may determine A(pc,Wi) numerically for given actions x by 
averaging the gravitational potential from Neptune over 7/2 at fixed W\, that is, 

A(x,w 1 ) = -^F{x,w 1 ), (23) 
a 

where 

q /-67T / \ r ■ r \ 

F(x, wi) — — I dy 2 \ - - - — — + constant; (24) 

67r J \ r n — r \T n r I 

\ 1 " 1 I'M / X,U>1 ,J/2 

the unimportant constant is chosen so that A = for circular orbits, i.e. F(xi,0,wi) = 0. 
It can be shown analytically that the contribution from the second (indirect) term in the 
integrand vanishes. 

We shall also write 

d 2 A M n ^ , . „ , , 2 d 2 F(x,w l ) , , 

A 22 = ^ = -—F 22 (x, Wl ) where F„{x,w 1 )=x 1 ^ ■ (25) 

Figure [l] plots the contours of F(x, 7f 1) at the resonant semimajor axis x± r , as obtained 
from equation fl24l) . The potential is singular for collision orbits, which for small eccentricity 
satisfy 

cos 7^1 = -. (26) 

ea 

The conditions (§) for stable zero-amplitude libration are satisfied if and only if w\ = w\ r = 71 
or Wi = and e > 1 — a n /a = 0.237. Figure |^ and |3| show examples of these two solutions, 
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plotted in a reference frame corotating with Neptune. Orbits of the first kind are similar to 
Pluto's, although with smaller libration amplitude (compare Fig. 4 of [Malhotra fc Williams 



|1997[ ). Orbits of the second kind (Fig. ^) were discussed by Malhotra (1996), who calls them 
"perihelion librators". We shall not discuss these further, since they do not appear to form 
naturally during resonance capture of initially circular orbits; moreover for e <i 0.35 they are 
likely to be unstable, since they cross Uranus's orbit and thus are subject to close encounters 
and collisions with that planet. 

Figure f| plots F(xi r , x 2 , 7r) and F 22 (xi r , x 2 , 7r); at the eccentricity of Pluto, correspond- 
ing to x 2 / (GM Q a) l l 2 = 0.0325, we have F(xi r , x 2 ,tc) = —0.313 and F 22 (xi r , x 2 , 7r) = 79.3. 
Thus, for example, the eccentricity oscillation in the separatrix orbit that marks the bound- 
ary between tadpole and horseshoe orbits has amplitude Ae = 0.007 (eq. |20[) and the period 
of libration of small tadpole orbits is 2ir/u = 9.1 x 10 7 yr (eq. |T9"D . 

For our purposes it is sufficient to work with the following numerical approximation to 
the resonant potential: 

~, x 0.584 + 0.130c, . , , 

(X ' Wl) = 1 + i 709e ln I 1 - 4.222 e cosH , (27) 

where x 2 = Xi r [l — (1 — e 2 ) 1 / 2 ]. This approximation is chosen to match the resonant po- 
tential at the resonant semimajor axis X\ = Xi r ; the dependence on the relative semimajor 
axes of Neptune and the Plutino is suppressed since the effects of this potential are only 
important near resonance. The logarithmic factor is chosen to reproduce the singularity in 
the resonant potential near the collision orbits defined approximately by equation fl26|) . The 
approximation formula also matches the analytic formula fl2l~D to 0(e 2 ) at w\ = it. 



Figure || shows the contour plot analogous to Figure ||] for the approximate resonant 
potential F, and the triangles in Figure § show F and F 22 . The agreement is very good, 
especially considering that errors are amplified by taking the two derivatives required to 
generate F 22 . 



3. Numerical experiments 

We follow the orbital evolution of Pluto and a Plutino in a simplified version of the Sun- 
Neptune-Pluto-Plutino four-body system that isolates the resonant potential from Neptune. 
Neptune is assumed to have a circular orbit that migrates outward according to the rule 
( IMalhotra 1993|) 

a n (t) = cif — Aaexp(— t/r), (28) 
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where cif = 30.17AU is Neptune's present semimajor axis, Aa = 6AU, and r = 1.5 Myr. 
Thus Neptune's initial semimajor axis is 24.17AU and the initial location of the 3:2 orbital 
resonance is 31.67 AU. 

The initial eccentricity of Pluto is taken to be zero and its initial semimajor axis is 
33 AU, as required so that its present eccentricity matches the observed value (eq. We 
followed 160 test particles, with initial semimajor axes distributed uniformly in the range 
[31 AU, 39 AU] and eccentricities distributed uniformly in the range [0, 0.03]. The inclinations 
of Pluto and the test particles are chosen randomly in the range [0, 3°] and their angular 
elements are chosen randomly from [0, 2ir]. Pluto and the test particles feel the resonant 
potential from Neptune, as defined by equations ( |2"B| ) and Q2"7|), but no other Neptune forces. 
The effects of the resonant Neptune potential on the orbital elements of Pluto and the test 
particles are followed using Lagrange's equations. 

The test particles do not influence Pluto or one another. However, they are subject to 
the gravitational potential from Pluto, 

B(x, y, t) = —GM P (— L- - r -^) ; (29) 
\\ L 1 p\ \ l p\ / 

the effects of this potential on the orbital elements of the test particles are followed using 
Gauss's equations. 

The evolution of Pluto and the test particles is followed for 0.45 Gyr or 10% of the age 
of the solar system. 

4. Results 

Of the 160 test particles, all but 12 are captured into the 3:2 resonance with Neptune, 
in the sense that their final semimajor axes are close to (3/2) 2,/3 a n and their eccentricities 
are near the prediction of equation ([5]), as shown in Figure |6]. The 12 particles that are not 
captured lie inside the location of Neptune's 3:2 resonance at the start of the calculation, 
(3/2) 2 / 3 x 24.17AU = 31.67AU, and would presumably be captured into other resonances 
if we used the full Neptune potential to work out their motion. We have verified this 
presumption by conducting shorter integrations (1 x 10 7 yr) using the same initial conditions 
but the complete Neptune potential. In this case all but 15 of the 160 particles were captured 
into the 3:2 resonance; the remainder were captured into the 4:3, 5:3 or 7:5 resonances. 

The behavior of the test particles in the 3:2 resonance (henceforth Plutinos) falls into 
the following broad classes: 
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Tadpole orbits (5 particles): these have longitude difference y\ — y p \ and differences 
in longitude of perihelion y p2 — y 2 that librate around the leading or trailing Lagrange 
point of Pluto (Figure [7|). (Note that the libration center for the orbit in this Figure 
is y\ — y p i ~ 100°, not 60° as implied by the analysis in §2.1. This discrepancy arises 
because Pluto has a high eccentricity, while our analysis is only valid for near-circular 
orbits. Similarly, the perihelion difference librates around w — w p ~ 300°, three times 
the difference in mean longitude as required by eq. 0.) The tadpoles show no evidence 
of chaotic behavior or secular evolution over the length of our integration. The analysis 
in §2.2 suggests that the maximum eccentricity difference for these orbits is Ae ~ 0.007; 
this requires in turn that their initial semimajor axes must have been close to Pluto's, 
as is seen to be the case in Figure []. 

Horseshoe orbits (19 particles): the longitude difference oscillates around 180°, with 
jumps in the Plutino eccentricity at the extrema of the longitude oscillation, as pre- 
dicted by the analysis of §2.1 (Figures ||, |]). The motion appears stable over the 
length of our integration although there are significant variations in semimajor axis 
oscillations during the course of the integration, and some horseshoes may evolve into 
transitional orbits over longer time intervals. 

Transitional orbits (43 particles): these show irregular behavior or transitions between 
libration and circulation of the longitude difference (Figure [TOD . When the longitude 
difference circulates, the particles are no longer protected from close encounters with 
Pluto. However, the particles remain in the 3:2 Neptune resonance in the sense that 
the resonant angle w\ continues to librate. 

Doubly transitional orbits (2 particles): Like transitional orbits, these show libration- 
circulation transitions in the longitude difference, but in addition they show irregular 
behavior in the resonant angle w±, leading eventually to a transition of w± from libration 
to circulation (Figure [11]). Although only 2 particles in our simulation exhibit this 
behavior, a number of others show growing amplitude in the W\ libration and will 
probably move into this class in less than the age of the solar system. Such orbits are 
normally short-lived since once W\ circulates, they are no longer protected from close 
encounters with Neptune. 

Irregular circulating orbits (17 particles): the longitude difference circulates through- 
out the integration. Pluto induces irregular behavior (Figure |12|), but the Neptune 



resonance is preserved in the sense that W\ continues to librate, at least over the span 
of our integration. 

Regular circulating orbits (62 particles): the longitude difference circulates throughout 



the integration, but the orbits appear fairly regular (Figure [13]). Generally, the orbits 
with larger eccentricity differences are more regular, because the encounter velocity 
with Pluto is higher so the perturbations from close encounters are smaller. 



These classes represent a sequence in eccentricity difference: the typical eccentricity 
difference |e — e p \ is smallest for tadpoles and largest for orbits unaffected by Pluto. Because 
the Plutino eccentricity is determined by the semimajor axis at the time of resonant capture 
(eq. HP, the classes also reflect the initial semimajor axes of the Plutinos: the tadpoles and 
horseshoes all have initial semimajor axes in the range 32.2 AU-34.2 AU (i.e. close to Pluto's 
initial semimajor axis of 33 AU). The transitional and irregular circulating Plutinos mostly 
have initial semimajor axes in the range 31.7 AU-36 AU, and the regular circulating Plutinos 
have initial semimajor axes concentrated in the range 35AU-39AU. 



5. Discussion 

Test particles captured into the 3:2 Neptune resonance (Plutinos) have a complex range 
of dynamical interactions with Pluto. The strength of the interaction depends on the dif- 
ference in eccentricity between the test particle and Pluto, and thus on the difference in 
initial semi-major axis if the initial orbits were circular and capture occurred through out- 
ward migration of Neptune. Plutinos are stable only if the eccentricity difference Ae is small 
(Ae £ 0.02 from Fi gure^), in which case the Plutinos librate on tadpole or horseshoe orbits; 
or if the eccentricity difference is large (Ae <> 0.06), in which case the longitude difference 
circulates but relative velocity at encounter is high enough that Pluto has little effect. Unsta- 
ble orbits at intermediate Ae can be driven out of the 3:2 Neptune resonance by interactions 
with Pluto, and thereafter are short-lived because of close encounters with Neptune. Thus 
we expect that the population of Plutinos has decayed over time, although determining the 
survival fraction will require integrations over the lifetime of the solar system using the full 
Neptune potential. 

The long-term behavior of orbits in the 3:2 Neptune resonance is central to the origin of 
Jupiter-family comets. The usual explanation is that slow chaotic diffusion and collisional 
kicks drive Plutinos out of the 3:2 resonance, after which they are subjected to close en- 
counters with the giant planets and eventually evolve into Jupiter-family comets ( [MorbideiT 



1997|) . Our results suggest that Pluto- induced evolution of Plutinos onto Neptune- crossing 



orbits may contribute to or even dominate the flux of Jupiter-family comets. 

Our results also enhance the motivation to obtain accurate orbits for Kuiper-belt ob- 
jects, and fuel speculation that the formation of the Pluto-Charon binary may be linked to 
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interactions between Pluto and Plutinos. 

This research was supported in part by NASA Grant NAG5-7310. We thank Matt 
Holman, Renu Malhotra, and Fathi Namouni for discussions and advice. 
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Fig. 1. — A contour plot of the resonant gravitational potential F(x,Wi) at the resonant 
semimajor axis, x% = X\ r (eq. |6]), as obtained from equation (pl[) . The ordinate is eccentricity, 
which is a proxy for x 2 = xi r [l — (1 — e 2 ) 1//2 ]. The contour levels are uniformly spaced at 
intervals of 0.01; positive contours are solid and negative contours are dotted. 

Fig. 2. — The orbit of a Plutino for 20,000 years in a reference frame corotating with Neptune. 
Initially the resonant angle w\ = 180° and the eccentricity e = 0.39. The high eccentricity 
was chosen to illustrate the character of the orbit; in a more realistic integration this orbit 
would be unstable because it crosses the orbit of Uranus. 

Fig. 3. — The orbit of a Plutino for 20,000 years in a reference frame corotating with Neptune. 
Initially the resonant angle w\ = 0° and e = 0.49. In a more realistic integration this orbit 
would be unstable, since it crosses the orbit of Uranus. 

Fig. 4. — The solid curves plot F(x lr , x 2 , vr) and (x 2 /x lr ) 3 ^ 2 F 22 (x lr , x 2 , it) as determined by 
the numerical integral ( ^4"] ) and equation (the multiplicative factor in front of F 22 is used 
because F 22 oc x^ 2 as x 2 — >• 0). The solid triangles and dashed line show the approximations 
F and F 22 given by equation ([27|). For reference, the eccentricity of Pluto, e = 0.253, is 
marked by an open circle at x 2 /(GM & a) 1 ' 2 = 0.0325. 

Fig. 5. — The contour graph of the approximation F (eq. ^) to the resonant gravitational 
potential. The contour levels are the same as in Figure |I[ 

Fig. 6. — The final eccentricity of the test particles as a function of their initial semimajor 
axis. The solid curve is the prediction of equation @, and the crossing point of the dashed 
lines shows the location of Pluto, computed with the same equation. The symbols denote the 
orbit classes discussed in §4: tadpole orbits (open circles), horseshoe orbits (solid circles), 
transitional orbits (crosses), doubly transitional orbits (open squares), irregular circulating 
orbits (starred squares), regular circulating orbits (solid triangles), orbits not captured into 
the 3:2 resonance (open triangles). 

Fig. 7. — A tadpole orbit. In the first panel "dpp" denotes the distance between Pluto 
and the Plutino in AU. Subsequent panels show longitude of perihelion, mean longitude, 
semimajor axis, eccentricity, orbital period ratio to Neptune, and the resonant angle of the 
Plutino; the detailed definitions of the orbital elements are in the text. 

Fig. 8. — A horseshoe orbit. 

Fig. 9. — A second horseshoe orbit. The spikes in A — A p arise because small short-period 
oscillations occasionally carry this angle past or 2tt. 
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Fig. 10. — A transitional orbit. 

Fig. 11. — A doubly transitional orbit. 

Fig. 12. — An irregular circulating orbit. 

Fig. 13. — A regular circulating orbit. 
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This figure "fig5.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9904424v 1 
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